NLO Polarized Chargino pair Production in Electron Positron Annihilation 
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We calculate the complete one-loop quantum corrections to the helicity eigenstate chargino pair 
production cross sections in polarized electron positron collisions, within the Minimal Supersymmet- 
ric Standard Model. We calculate the non-QED corrections using the helicity amplitudes formalism, 
and Dimensional Regularization to deal with ultraviolet divergences. We calculate QED corrections 
using the dipole subtraction formalism to extract soft and collinear divergences in Bremsstrahlung, 
canceling them with the infrared divergences from virtual QED corrections. We show numerical 
results for the Focus Point scenario in mSUGRA, where we find important quantum corrections for 
differential cross sections with definite chargino helicities. 
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I. INTRODUCTION 

In supersymmetry the fermionic partners of charged Higgs and W gauge bosons, the higgsinos and 
winos, mix to form a pair of charged fermions called charginos, \f liD- I n many supersymmetric scenarios 
they are light enough to be produced at the LHC and ILC, although they have not yet been observed (2J, 
6, 7]. The discovery of heavy charged fermions (8|, |9|] not coupled through strong interactions 
is not proof of supersymmetry, though. Precise measurements of their masses and couplings Jiol should 
exhibit the supersymmetric prediction that SM couplings have their mirror coupling with superpartners. For 
example, WW Z coupling should be equal to the Z coupling to a pair of winos. This, together with the 
experimental precision expected at the ILC, leads to the necessity to have precise theoretical calculations 
for the observables in order to properly compare them with the experimental results 
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positron collisions |15 



and chargino decays (23 
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14H . chargino production cross section in electron 



chargino production cross section in hadron colliders [22], 



2411 . are well documented. Here we improve the collective knowledge by pre- 



senting a complete one-loop calculation to the production cross section of polarized charginos in collisions 



between polarized electrons and positrons 125 
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28]. We stress the fact that this is the first time com- 



plete NLO corrections to the production of polarized charginos are calculated. These quantum corrections 
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are an effect of the rich electroweak properties of the MSSM, and may provide a window to the Parity 
violating structure of this supersymmetric model. The full potential of this calculation will become avail- 
able when quantum corrections to polarized chargino decays are incorporated, as indicated by the tree-level 



analysis of the importance of spin correlation between chargino production and decay H29TI . 

In order to handle the full one-loop calculation we use form- factors to organize the bubble and triangular 
diag rams, and Q-charges to merge form-factors with box diagrams. We use the helicity amplitude formalism 



[25, 



26, 



2711 to keep track of the electron and positron polarizations and the chargino helicities. To calculate 
the QED corrections and cancel the Infrared (IR) divergences we use the dipole subtraction formalism 
introduced in [30]. In this method the IR divergences are isolated analytically through the introduction 
of auxiliary subtraction functions. In this way, an IR cutoff is not needed, avoiding unpleasant numerical 
problems that emerge when a cutoff is used to regularize the infinities. At the end, the three body phase 
space integral is performed using a Monte Carlo method. 



II. CHARGINO PRODUCTION AT TREE LEVEL 



In the MSSM, the scattering e + e — > x + X arises at tree level through 7 and Z in the s-channel, and 
through v e in the t-channel, as indicated in Fig. Q] In this article we include the effects of the polarization 





FIG. 1: Three level diagrams for chargino pair productions at the ILC. 



of the electron and positron, and the chargino helicities. The notation for the four-momenta is as follows, 



e + (p2)+e (pi) -» xt (fa) +X a (h) 



(1) 



where the convention is that pi and P2 are incoming and k\ and k 2 are outgoing four-momenta. The Z 
coupling to two electrons is, 

9 



y Zee 



2cw 



7 M [9v ~ 9 Alt] 



(2) 



where c\y = cos 9w, % = sm0\y, and 9w the weak mixing angle, with sin 0% S (m z ) = 0.23122 feUl 
Also g v = —1/2 + 2s^ and g A = —1/2 are the vector-axial couplings. The SU(2) gauge coupling 



constant g is related to the positron electric charge by g = s\ye, with a{mz) = e 2 /Air = 1/127 



We 



also work with m z = 91.1876 ± 0.0021 GeV and T z = 2.4952 ± 0.0023 GeV |33] for the mass and width 
of the Z boson respectively. The Z-coupling to a pair of charginos is, 



-,ab vfi 
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n r U ab ■ 



(3) 
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where we use the notation in ref. [ 1] for the 0' L and 0' R couplings. 

The photon couplings to a pair of electrons and a pair of charginos is simply given by their electric 
charge, 



->ab u,0 

y 7XX 



ab 



(4) 



Note that our convention is that a positive chargino is the particle and a negative chargino is the anti-particle. 
The couplings between electrons charginos and sneutrinos are, 

1 



->&+,0 



1+75^ 

-gVbi — n — c > 



75 



(5) 



where the matrix V is one of the matrices that diagonalize the chargino mass matrix, in the notation of |l|], 
and we are assuming CP is conserved and couplings are real. The matrix C is the charge conjugation matrix, 
and appears due to the clashing arrows in the sneutrino Feynman diagram in Fig.Q] Only experimental lower 
bounds for the sneutrino mass have been set, of which we mention > 94 GeV Q3[] . 



III. QUANTUM CORRECTIONS 

We regularize divergent diagrams using dimensional reduction DR. In each graph, divergences are 
contained in the parameter 

2 

A = - + In47r-7E (6) 

4 — n 

where 7^ is the Euler-Mascheroni constant, and n is the number of space-time dimensions, which is defined 
to be n = 4 — 2e, with the limit e — > to be taken at the end of the calculation. We call the renormalization 
subtraction point Q, and it is such that the parameters that define the tree-level cross section are promoted 
to DR running parameters, evaluated at the scale Q. We perform our calculations from first principles, 
not relying on loop calculating packages, as a way to find an independent result from groups that use these 
packages. 

We have two types of divergences. The ultra-violet (UV) divergences appear at large internal momenta 
running in the loops, and they are to be cancelled with counterterms introduced by the renormalization 
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procedure. We also have diagrams with infrared (IR) divergences, which appear at low internal momenta. 
These IR divergences are cancelled by bremsstrahlung, i.e., real photon emissions when the photon is soft 
or collinear. 

We organize the calculation combining two formalisms, the form-factors and the helicity Q-charges. The 
form-factor formalism is specially useful to organize bubbles and triangle graphs, introduced in ref. I15I in 
the context of the quark-squark approximation, in which only quark and squark loop corrections are taken 
into account. The helicity amplitude formalism is necessary when boxes are included, and was introduced 
in ref. 12511 where non-QED boxes were considered. It is specially designed to keep track of the polarization 
of the electrons and positrons and the helicity of the charginos, in the context of full one-loop calculation 
we are dealing with in this paper. We describe and combine both formalisms below. 

A. Form-Factor Formalism 

n 

In ref. IU5I1 we calculated the radiative corrections to chargino production in electron-positron collisions 
in the approximation where only loops from quark and squarks were considered. A form-factor formalism 
was developed, and here we generalize it to include triangular graphs in the gauge boson vertices to electron- 
positron pairs as well. 

We organize the triangle graphs with the help of form-factors that define Green's functions where the 
fermionic lines are on-shell, but the bosonic line are not. In the case of Z couplings to a pair of charginos, 
they are given by 

Sz xx - F z xx T g + Zxx 7 2 Zxx ~/i 2 Zxx T/i 2 ( } 

and similarly for G^ X x replacing Z — > 7 in the form-factors. The v and s labels refer to vector and scalar 
couplings respectively. We define fc_ = k± — ki, which is normalized by ^fs in order to make the form- 
factors dimensionless. Note that in ref. iQ] we used a slightly different but easily related definition for the 
form-factors. 

In the case of e ± -sneutrino-chargino vertices the form-factors are, 

fe+ = + (l+7 5 ) r C a- = p - r -l (1 ~ 7 5 ) fg , 

y i/ex vex 2 ' x uex 2 

where C is the charged conjugation matrix. In eq. ©, the ± notation refers to the positive (negative) 
outgoing chargino. 

Here we extend the formalism in ref. I15I to include form-factors in the electron-positron vertex to a 



gauge boson, 



/-•H IPVR „,U (1 7 ) . rpvL _.It(l 7 ) /Q\ 

G Z ee = F Zee T o h ^ Zee T ^ (9) 
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and similarly for the photon. Only vector form-factors are needed because scalar ones disappear when 
electrons and positrons are taken on-shell. Also note that in the quark-squark approximation used in ref. Jl5l1 
none of the Zee form-factors receive quantum corrections. 

We list below the form-factors that do not vanish at tree-level, 

rpvR,0 _ _9_r\lR t?vL,0 _ 9 ^/L 

r Z XX ~ c w U ab > r Z X X ~ c W b 

pvRfl _ r rwLfi = r 

= -gV u , F-^l = gV al (10) 

Fzef = -^(f V ~ 91) , Fief = + 9%) 

t?vR,0 _ r?vL,0 _ 

"fee ° 3 7ee 

We highlight a few details about these tree-level form-factors. First, the difference in sign between F 7XX 
and F 7ee is due to the fact we consider electrons and positive charginos as particles, and positrons and 
negative charginos as anti-particles, and of course F 1XX is proportional to 5 a b due to electromagnetic gauge 
invariance. Also, we note that the form-factors F„ ex indicate that the wino component of the chargino is 
the one that couples to the sneutrino, diminishing considerably the t-channel contribution in the case of 
higgsino like charginos. In contrast, since, 

0$ = -U al U bl - \U a2 U b2 + s 2 w 8 ab , 0% = -V al V bl - \V a2 V m + s 2 w 5 ab , (11) 

we see that both the wino and higgsino components of the charginos couple to the Z gauge boson. Finally, 
we note that since sin 2 6y/ is approximately equal to 1/4, the Z vector coupling gy to a pair of electrons 
is much smaller that the axial one g\. As a consequence, F^f F^f are nearly equal in magnitude but 
with opposite sign, which in turn implies a cancellation in the right handed electron cross sections. This 
cancellation will become more clear in the following section. 

B. Helicity Amplitudes Formalism 

The helicity amplitude formalism is a very systematic and economic way to organize virtual corrections 



through the so called Q-charges, introduced in 12511 . Here we calculate also QED corrections not included 
in the above reference. 

Consider the scattering amplitude for electron-positron into a pair of charginos. On the one hand, the 
electron has a polarization a = R, L and momentum p\, while the positron has opposite polarization and 
momentum p 2 . On the other hand, the negatively charged anti-chargino has a mass m Xa , momentum k%, 
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and helicity Ai, while the positively charged chargino has a mass m Xb , a momentum k 2 , and helicity A2. 
The dimensionless scattering amplitudes are written as, 

A% M = - s L lxa Q» i (k 2 ,\ 2 \r i \k 1 ,\ 1 ). (12) 

where there is no sum on a. The normalization is such that the differential scattering cross section is 



d cos 32 7r s 



and X(x, y, z) = x +y + z — 2xy — 2xz — 2yz. The leptonic matrix elements are, 

L m = v(p 2 h ^- 2 M n(pi) (14) 

with the +(— ) sign for a = R(L). Note that there are no more structures because we are considering 
massless leptons. The chargino matrix elements are calculated with the help of five Dirac structures P, i = 

1>2 = (1±7 5 ) 
2 

r s,4 = yil±2f) (15) 

T 5 = -ia vp 

Note that these structures can have undisplayed Lorentz indices which are in turn contracted with cor- 
respondingly undisplayed Lorentz indices in the dimensionless tensor coefficients Q^. The contribution 
from any Feynman graph to such an amplitude can always be expressed in this form by making a suit- 
able Fierz transformation where necessary. The tensor coefficients can be further reduced to scalar 
Q-charges Qf, 



= 


Qf ktl^l 






= 








QZ = 




kl 


P-/s 


V Q 4 — 


Q%^ v +Q2,2 


tit 


P-/s 


n^f — 









(16) 



where again, the scalar Q-charges Qf may have an extra numeric label, and they are dimensionless. In 
the above equations we have defined kt = (fef — k% ) and p^_ = (pi — p 2 )- Any other structure can be 
expressed in terms of the above quantities, by exploiting the fact that the leptonic current is conserved and 
that the matrix elements of P are taken between on-shell chargino states. 
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In this article we concentrate on the case of chargino pair production of equal mass. We define k 



e m 



| A;i | = | k 2 | the magnitude of the 3-momentum of each chargino in the CM-frame, 9 the scattering ang 
the same frame, and v = Ikj^fs its velocity. Helicity amplitudes are given for the general case in ref. [25], 
and here we list them for the case m Xa = m Xb . Helicity amplitudes are denoted by A® 2 ^ , where a = L,R 
is the polarization of the electron, and A2A1 are the helicities of the chargino and anti-chargino respectively. 
For left handed electrons we have: 

A+ + = -Qiv(l-v)sm8 + Q%v(l + v)sm9 



+(Qh + 24i) V 1 -v 2 sm6 - (Qg 2 + Q^v^l-v 2 sin 9 cos 9 

-2Q^v sin0 + 4Qf 2 sin0 (17) 



A%_ = -Qfi(l + u)(l + cos 9) - Q% 2 v{l + v) sin 2 9 
-Q4i(l-«) (1 + cos0) - S£ 2 «(l - v) sin 2 9 

-4Q% 2 Vl-v 2 {l + cos9) (18) 

At+ = +Sfi(l - v) (1 - cos 9) - Q% 2 v(l-v) sin 2 9 
+Qfi(l + v) (1 - cos 9) - Q£ 2 v (1 + v) sin 2 9 



+4Qf 2 V1^2(l-cos0) (19) 

= -Qf't;(l + t>)siii0 + Q£u(l-u)sin0 

-(Q31 + Q4i)Vl sin0 + (Qf 2 + Q^^l-w 2 sin cos 9 

-2Q^v sin 9 - AQ^ 2 sin 9 (20) 

where v is the chargino velocity given by 

/ 4m?. 

(21) 
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The helicity amplitudes for right handed electrons are: 

A++ = -Qf v(l -v)sm9 + Q$v(l + v)sim 



+(Qfi + Qfi)Vl~v 2 sm6 - (Qf 2 + QlD^V 7 !-^ 2 sin^ cos 

+2Q^ u sin - 4Qf 2 sin 9 (22) 

.4+_ = +Qfi(l + u) (1 - cos#) - Qf 2 v(l + v) sin 2 
+Qfl(l - v) (1 - cos 0) - Qf 2 v (1 - v) sin 2 



-4Qf 2 >/l -v 2 (1 - cos 0) (23) 

A R + = -G&(1 - u) (1 + cos#) - Qf 2 u (1 - u) sin 2 9 
-Qfl(l + u) (1 + costf) - Qf 2 v(l + v) sin 2 

+4Qf 2 \A - v 2 (1 + cos 9) (24) 
„4 R _ = -Qf v(l + v)sm9 + Q^v(l-v)sin9 



-(Qfi + Cfi)-s/l -^ 2 sin0 + (Qf 2 + Qf 2 )u Vl-v 2 sin^ cos 

+2Qfj u sin + 4Qf 2 sin (25) 

The simplicity of these expressions is striking. All the quantum corrections are concentrated into the Q- 
charges. To find differential cross sections for specific chargino helicities, we just have to square the cor- 
responding amplitude above, and feed it into eq. (TT3T ). If we sum the squared amplitudes over chargino 
helicities we obtain from the expressions above the polarized differential cross sections when the beam is 
100% polarized. But within this approximation, we do not include spin correlation. The full potential of the 
result of this article will be realized once the decay chain is added, calculating the chargino spin conelation 
between production and decay [29], and its effect on observables like angular distributions of final decay 
products. 

It is instructive to analyze the tree-level approximation with the help of the Q-charges. The tree-level 
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expressions for the Q-charges are, 

e 2 d ab 

e 2 S ab - g 2 V b i S 2 Ki 
t — m~ 

e 2 S ab (26) 
e 2 5 ab 

This coincides with the expressions given in ref. [10] (after allowing for the fact that in ref. [10] the nega- 
tively charged chargino is taken to be the particle and the positively charged one the antiparticle, whereas 
our convention is vice versa). Note that the only dependence on the scattering angle 9 of the Q-charges at 
tree level comes from the i-Mandelstam variable from the sneutrino contribution, which is always negative. 

With the help of eq. d26l ) it is easy to understand why the right handed electron cross sections are smaller 
than the left handed ones. Since g v is very small and g e A is negative, the coefficients for the O' couplings are 
positive for left handed electrons and negative for right handed ones. On the other hand, for equal charginos 
the O' couplings reduce to, 

O'fi = | sin 2 <j> R - c 2 w , 0[\ = | sin 2 </> L - c 2 w (27) 

where 4>r is the rotation angle that defines the rotation matrix U, and analogously for the angle ob L and the 
matrix V. In this way, the 0' aa are always negative, implying that the photon and Z contributions interfere 
destructively for right handed electrons and constructively for left handed ones. 



L '° 

^31 



o L,0 



Q R,0 
^31 

R,0 
41 



Q 



In 2 



2c 2 



2r 2 



+ 9%)- 



m 



n' R 

2 U ab 



+ 9 e A> 



-O 



7717 



iL 

ab 



)V - 9a) 



-O 



iR 

ab 



2r 2 



(g e v-9A> 



-o 



rrir 



IL 

ab 



IV. NON-QED VIRTUAL CORRECTIONS 

In this article we present separately QED corrections (bremsstrahlung and virtual photonic graphs), from 
the rest of the virtual graphs. The later are analyzed in this section. We distinguish three kinds of diagrams 
which contribute to the total cross section: bubble (two point functions), triangular (three point functions), 
and box (four point functions) diagrams. We work in the Feynman gauge £ = 1. 

Bubble diagrams can be easily inserted into the form-factors defined before. We consider first the two- 
point Green functions with gauge bosons in the external legs. The two gauge bosons, which we call V and 
V' , are off-shell, and in our case they correspond either to neutral Z or 7. In Fig. [2] we see the diagram for 
this Green function, whose Lorentz structure is defined with the help of two scalar functions A and B, 

^ vl {p 2 ) = t[A vv ,{p 2 )g^ + B vv ,(p 2 )p^], (28) 



10 



FIG. 2: Vector boson two-point Green function. 



Complete prototype for bubble diagrams will be given elsewhere, while the A function can be found in B25I1 



The photon self-energy vanishes at zero momentum by virtue of gauge invariance so that after subtrac- 



tion in the DR scheme contributes to the photon form-factor according to 

Ai^f ) = -e^±5 ab (29) 
where a and b refer to the two species of charginos produced. The photon- Z mixing is also subtracted in 



the DR scheme and contributes to the Z form-factor 



and to photon form-factors: 



AFfJ^ = -e^±5 ab (30) 



1XX cw s- to| 

The Z-boson self-energy is regularized with a subtraction at s = m 2 z : 

AF ^X =—°a b J^T W 

Another important two-point Green function is the sneutrino self energy, whose diagram is in Fig. [3] 
Since it is Lorentz scalar, it is represented by only one function A, 



FIG. 3: Sneutrino self energy. 



^ e (p 2 ) = i[A^ e (p 2 )]. (33) 

The contribution to the sneutrino form-factors is obtained after regularizing with a subtraction at t = to~, 
with the following result, 



AF ,e X ~ 2^(»)1 ^1 (34) 



11 




FIG. 4: Chargino two-point Green function. 



This guarantees that the parameters mz and rriy respectively refer to the physical (pole-)masses. 

The chargino two-point function can be seen in Fig. [4] where the external legs correspond to any of the 
two charginos. Its Lorentz decomposition involves four scalar functions, 



AUp 2 ) + BtAp 2 )j> 



;i+7 5 ) 



(l-7s) 



(35) 



The chargino self-energy and mixing contribute to form-factors in a more complicated way. Since these 
are external particles we have insisted that the subtractions are performed on-shell, so that the renormalized 



25b. 



chargino fields are indeed physical fields. Details can be found in ref . |15 , 

Also embedded into form-factors are triangular diagrams, grouped into Zxx> Zee, 7XX> 7 ee an d e X& 
three-point functions. Form-factors are in turn easily incorporated into Q-charges. On the one hand, the left 
handed Q-charges that receive contributions from form-factors are, 



AQ31 
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2c w 
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2c w 

9 
2c w 

9 
2cw 



rrir 



(9v+ 91)^2 AFI L xx + eAF% x 



rrir 



(36) 



(9v+9a)- 



rrir 



(g e v + 9 e A> 
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U ab — 

C\V S 
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rrir 



-AFf ee + eAF^ x - e5 ab AF^ e 



•yee 
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AFr xx + —Of b - 
z c w s 



rrir 



-AF v z l + eAF^ xx - e5 ab AF? L 



■yee 



1 



to- " A 2" ^t — m- 



tAF- 



+ 



On the other hand, the right handed Q-charges receiving contributions from form-factors are given by, 
9 
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AQ[ 



AQ2 
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(g e v-9A> 
(g e v-g A > 
(g e v-g e A> 



-AFz xx + eAFf R 



rrir 



ixx 



m 



2 AF Z XX + eAF ' 



IXX 



(37) 



to: 



-AFt xx + — O 



—<1' R — 

Cw S 



rrir 



-AF vR e + eAK 



ixx 



e5 ab AF^ e 



m' 



A vp vL _i_ ^ n iL * 

' IXt Z X X + ~°ab 



C W 



rrir 



-AF vR p + eAF% x - eS ab AF^ R 



■ -fee 



Finally, box diagrams are directly incorporated into Q-charges with prototype diagrams in ref. [25], and 
more general diagrams that will be shown elsewhere. 
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As opposed to bubbles and some triangular diagrams, box diagrams are UV finite. Ultraviolet diver- 



gences that occur in a few of the triangle graphs are subtracted in the DR scheme with subtraction point Q. 
Therefore, apart for the masses which are taken to be physical and the weak-mixing angle whose renormal- 



ization is described above, all other parameters are to be considered to be in the DR at the scale Q in the 
MSSM theory. This means that the translation of the values used here to those directly extracted from exper- 
iment, such as neutral current neutrino scattering cross-sections or the measured fine-structure constant will 
be slightly different from that of the Standard Model (without the supersymmetric partners). For example, 
the treatment of the photon- Z propagator system, described above, guarantees that the propagators only 
have poles at zero and Mz, but there is still some remnant of photon-Z mixing at these poles. We have 
checked numerically that the effect of a further subtraction of the photon-Z mixing propagator to remove 
this mixing has a negligible numerical effect on our results. Furthermore, the input SUSY parameters cho- 
sen are assumed also to be the corresponding values renormalized in this scheme at the same scale. We 
expect the sensitivity to (reasonable) changes in renormalization scheme to be genuinely of order aw / ft 
and to have no significant effect on our numerical results. The numerical value we use for the subtraction 



scale in this article is Q = 1 TeV as suggested by the SPA Project [32]. 



V. QED CORRECTIONS 

In the previous section we discussed the non-QED virtual corrections to our chargino production process. 
In this section we discuss the treatment of the remaining QED corrections separately because special care 
must be taken with IR divergences present in bubbles, triangles, and boxes involving photons. These IR 
divergences appear in graphs with virtual photons at low momenta of internal particles in the loop. They 
cancel from the cross section with correspondin g IR divergences comming from real photon emission from 



the external fermions, i.e., bremsstrahlung B34l 



35, 



36], which appear in both soft and collinear photons. 
As a remanent, large corrections proportional to ln(s/mg) remain, known as Leading Logarithms (these 



are initial state radiation corrections) and common to any of such similar process. For bremsstrahlung 
corrections we do not use the Q-charges formalism, calculating directly the corrections to the amplitude 
squared, using the algebraic manipulator FORM [38]. 



A. Bremsstrahlung corrections 

Diagrams contributing to the process e + e~ — ► X + X~l w& depicted in Fig. [5l where the bubbles repre- 
sent the three different tree-level channels (7, Z, and u e ). 
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X T e 



X e 




X T el 



X e" 



FIG. 5: Real photon emission. 



The total cross section o rea i for a real photon emission can be written as follow, 



space 



(38) 



where ci<I>3 is the three-body phase space and A 7 is the photon polarization. We implement the dipole 
subtraction formalism introduced in [30] in order to analytically isolate the IR divergences coming from 
the emission of real photons, and cancel them from the also analytically extracted divergences from virtual 
photon diagrams. Once IR divergences are removed, the QED corrections are IR free and are calculated 
numerically. The dipole subtraction formalism conveniently avoids numerical problems arising when a 
cutoff is used to regularize the infinities. 

The dipole subtraction method entails the introduction of an ad hoc auxiliary function |.M su b| 2 which 
becomes equal to \M- re ai\ 2 in the soft and collinear limits, 



\M 



sub\ 



^2\M rec a\ 2 , as pi ■ k 



0. 



(39) 



where pi are the four-momenta of the massless external fermions and k is the four-moment of the photon. 
Adding and subtracting this function into eq. (|38T > we find, 

(40) 



(Treal = J d<S> 3 \M rea l\ 2 ~ \M sub \ 2 ) + J \M suh \ 



where the first integral is finite and is performed numerically without the serious instabilities found in cutoff 
methods. The second integral can be done after an analytical extraction of the IR divergences, regularized in 
I30I1 by introducing a non-zero photon mass m 7 , and an electron mass m e . Once the (universal) contribution 
proportional to ln(m e ) arising from initial-state radiation is removed, the soft and collinear divergent terms 
can be mapped into pole terms in dimensional regularization using the mappings 
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These IR divergences cancel against IR divergences from virtual diagrams (also treated by dimensional 
regularization). 

With the process described above, the IR divergences from bremsstrahlung diagrams are, 



aAi A2 
real, IR 



2vr 



2 5 2, (Q z \ 2l + v\ 1-v 4, mf - it 

+ - + - In — + In + - In — * 

- c £ Vc / 1 + v e m~ — t 



v 



\ 



(7. 



2^2 



(42) 



which as we see, are proportional to the tree-level chargino pair production cross section. In the above result, 
the double pole comes from simultaneous soft and collinear divergences, the simple pole not proportional 
to any logarithm comes from collinear divergences, while the simple pole proportional to logarithms comes 
from soft divergences. 

B. Virtual Corrections 

Virtual QED quantum corrections at NLO in a e are given by, 

<rt A2 = £ fd*2 2Re((M«tt TM^) (43) 
where * s tne iree level amplitude, ■M^ r 1 t X2 contains all QED virtual corrections, and the photon 



coupling to fermions is factor out in a e /(47r). In DR, the IR divergent contributions from virtual QED 
correction can be written as, 



a. 



a.\\ A2 
virt, IR 



'2vr 



2 5 2, Q 2 \ 21 + v\ l-v 4, m| - u 
+ _ + + l n + _ l n 

e z e e V s / 1 + u e — 1 

X. 



a. 



0A1 A2 
2^2 



(44) 



where details will be given elsewhere. Clearly from eqs. (1421 and (l44l . virtual and real IR divergences 
cancel each other. 



VI. NUMERICAL RESULTS 



As a working example we concentrate on the mSUGRA Focus Point known as SPS2 [39]. This scenario 
is characterized by, 

m = 1450 GeV, m 1/2 = 300 GeV, A = GeV, tan (3 = 10, fj, > 



and we use the code ISAJET for the running from GUT to weak scale II37I1 . The low energy soft parameters 
calculated this way are fed into our code. The integration over phase space is performed with a MonteCarlo 
technique. Some relevant low energy parameters and masses are shown in the following table, 
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parameters (GeV) 


masses (GeV) 


Mi = 121.9 
M 2 = 235.7 
H = 222.2 


m°, = 172.1 

xf 

m° + = 297.0 

xi 

m^o = 111.5 

Xi 

m/> e = 1454 



where the charginos, neutralino, and sneutrino masses are given at tree-level. We correct the chargino 
masses at one-loop, and in order to compare the corrected cross section without moving the threshold, we 
define a tree-level scenario where the values of Mi and \i are tuned such that this tree-level masses are 
equal to the one-loop corrected chargino masses. These new tree-level parameters and one-loop corrected 
chargino masses are in the following table, 



parameters (GeV) 


masses (GeV) 


M 2 ° = 242.9 
H° = 219.7 


m~+ = 173.7 
m~+ = 300.0 

X? 



We report below only on production cross sections of light charginos, with a mass of m + = 173.7 GeV as 

Xi 

shown on the table. 

As we mentioned before, the renormalization scale we use is Q = 1 TeV motivated by the SPA Project 
0211 . Electroweak observables sin 2 6^r S {m z ) and a(mz) mentioned in section II must be run up to the scale 
Q = 1 TeV We find sin 2 6^(0) = 0.2395 and a(Q) = 1/125, which translates into g(Q) = 0.6479 and 
g'(Q) = 0.3636, as reported in 132J]. 

We separate the corrections in the following way: by "MSSM" corrections we mean all loops that do 
not include photons, by "QED" corrections we mean all loops including photons (virtual QED) plus all 
bremsstrahlung diagrams, and by "££" we mean the leading logarithms, which are universal to all this kind 
of processes. 

Our definition for the "££" contribution is, 



a 



aAiA2 
U 



-M4 

7r V mi 



dx 
o 1 



l + x 2 



M aXlX2 



l 2~>2 



(xp) 



d$ 2 (i; 



(45) 



where the term In ^^-J is precisely the large logarithm that gives the name to the whole leading logarithm 
contribution. The function (1 + x 2 )/(l — x) is the splitting function, which is related to the probability 
of finding an electron or a positron with a momentum xp before the emission of a photon, which takes a 
momentum (1 — x)p. The differential d<&2{x) is the 2-body phase space with one incoming 4-momentum 
multiplied by x, which translates into s — > xs. 
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A. Unpolarized Total Cross Section 
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FIG. 6: Unpolarized total cross section as a function of ^/s for the SPS la scenario. In the upper frame we compare 
the tree-level (red-dashed), with full NLO (black-solid). In the lower frames, the different curves correspond to full 
NLO corrected (black solid), MSSM (blue dashes), QED (red dot-dash), and 11 (green solid). 



We start by comparing our unpolarized results with the ones reported in 11711 . In this reference the result 
for NLO corrections to the unpolarized cross section was reported for the scenario known as SPS la. In that 
case, complete NLO quantum corrections are very large and negative, strengthen by a large contribution 
from sneutrinos, due to a small sneutrino mass and large sneutrino-neutrino-electron coupling. Our predic- 
tions for the total NLO cross section for light charginos in SPS la, whose corrections vary from —27% at 
^fs = 400 GeV, to —7% at yfs = 1400 GeV, are in reasonable agreement with the ones reported in I17I1 
within ±2%, as can be seen in the lower-left frame of our Fig. [6] Similarly, the total NLO cross section 
for heavy charginos vary from —22% at yfs = 980 GeV to —8% at yfs = 1400 GeV, and they are also in 



agreement within ±2% with 11711 . 
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Of course, the magnitude of quantum corrections depends on the supersymmetric benchmark chosen. 
We already mentioned that our working scenario is benchmark SPS2, also known as Focus Point. In Fig. [7] 
we show the one-loop unpolarized total cross section as a function of the center of mass energy y/s. In 
the upper plot we have the total cross section in the Born approximation (red dash line), and the NLO 
corrected total cross section (black solid line), for light and heavy charginos. We see that in this scenario 
NLO corrections to the total cross section for light charginos, change sign at i/i as 1250 GeV, while for 
heavy charginos the correction is always negative in region shown in the graph. In the lower frames of 




FIG. 7: Unpolarized total cross section as a function of y/s. The different curves correspond to full NLO corrected 
(black solid), MSSM (blue dashes), QED (red dot-dash), and tl (green solid). 

Fig. [7] we show details of the NLO corrections, as in the previous figure. The magnitude of the corrections 
rise sharply as the energy approaches the threshold, although this percentage correction acts on a total cross 
section that approaches zero in this limit. 

We note that the corrections in the case of Focus point SPS2 considered in this paper, are considerable 
different from those of benchmark SPSla. The MSSM correction (excluding the pure QED corrections) are 
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consistently lower (negative) in the case of SPS2. The major difference comes form the 1.1. part of the QED 
corrections which, at high energies, are nearly three times as large as the case of SPSla. This substantial 
difference is in turn dominated by a difference in the initial radiation correction to the sneutrino exchange 
part of the amplitude. For SPS2 the sneutrino mass is much larger than for SPSla, so that its contribution 
to the tree-level differential cross-section is very insensitive to scattering angle, and negligible in front of 
s-channel contributions. On the other hand, sneutrino contribution in SPSla scenario is large, and interferes 
destructively with s-channel contributions. All in all this leads to a light chargino total correction for the 
SPS2 point which is positive for y^i > 1250 GeV, whereas for SPS la, the total correction remains negative. 

B. Polarized Differential Cross Sections 

In this section we show our results for one-loop corrected differential cross sections for the production 
of light charginos with definite helicity from polarized electron positron collisions, in the Focus Point SPS2 
scenario. We choose to show our results for the differential cross sections as a function of xf transverse 
momentum p? (rather than the scattering angle) defined as, 

Pt = ~(Pi • fe)(p2 • k 2 ) - ml+ (46) 

s *l 

where the center of mass energy is yfs = 1000 GeV. In the case there is no photon in the final state, this 
reduces to px = yfs v sin 9/2, where v and 9 are \t velocity and scattering angle. 

The differential cross section we are working with is related to the differential cross section as a function 
of cos 9 in the following way, 

da 2 tan 9 da 

1— = a (47) 

dpT V s v cos " 

which partly explains the tendency for the cross section da/dpx to grow with px, as we go from 9 = to 

9 = -k/2. 

In the following six figures we show quantum corrections to polarized chargino production cross section 
as a function of transverse momentum. All of these figures have the same structure and differ in the chargino 
helicities and the electron polarization. In all plots we show px > 100 GeV because it is the most interesting 
region, with larger cross sections and better chances of differentiation from background. 

In Fig. [8] we have da\~ /dpx, ie., for a chargino with positive helicity, an antichargino with negative 
helicity, and a left handed electron (right handed positron). On the left we have the differential cross 
section as a function of pr, where we compare the born approximation, the one-loop corrected cross section 
including MSSM and QED corrections, and the corrected cross section including MSSM, QED, and leading 
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du L + /dp T in Focus Point Scenario 




Pt (GeV) p T (GeV) 



FIG. 8: da+-/dp T as a function of px in the Focus Point scenario within mSUGRA. 

logarithms. QED+MSSM corrections are moderate and negative, varying between -1% and -18%, as can 
be seen also in the lower right frame. Leading logarithms are very large, and increase the differential cross 
section even more, except at large pr, where il are negative. The largest positive corrections are obtained for 
low transverse momentum, px ~ 100 GeV, while the largest negative corrections appear at ~ 440 GeV. 
This choice for the chargino helicities and electron polarization gives the highest values for the differential 
cross section, varying between 0.2 and 2.1 fb/GeV. In the upper right frame we show the MSSM and QED 
corrections separately from each other, and compared with the Born approximation. We see that the QED 
corrections (not including leading logarithms) are smaller than 1%. 

In Fig. [9] we have da^ + /dpT, which corresponds to the case with both chargino helicities reversed with 
respect to the previous case. The differential cross section is comparable, varying between 0.1 and 1.4 
fb/GeV, as we see in the left frame. In the lower right frame we observe that QED+MSSM corrections are 
larger in magnitude than the previous case, and also negative. Note that here the relative magnitude of the 
corrections decrease with px, as opposed to the previous case, where they increase. Another difference with 
the previous case is that QED corrections are larger, but smaller in magnitude than MSSM corrections, as 
seen in the upper-right frame. This is a case where QED quantum corrections should not be neglected. The 
Leading Logarithms are usually positive, but become negative at large px, owing to the fact that at large 
Pt, x is required to be larger than 4p T /s, thereby reducing the x-dependent phase-space factor &2{x) in 
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da L + /dp T in Focus Point Scenario 




p T (GeV) p T (GeV) 

FIG. 9: d(T^ + /dpx as a function of px in the Focus Point scenario within mSUGRA. 

eq.(l45T). All the above differences between the corrections to o\~~ and are interesting manifestations 
of the parity violation properties of the MSSM. 

In Fig. [10] we see the equal chargino helicity production cross section. CP invariance, which we assume, 
enforces the equality of da^ + /dpx and da^~ /dpx, and it serves as a check of our calculations. This cross 
section is much smaller, varying from less than 0.001 fb/GeV at low px up to almost 0.2 fb/GeV at large 
PT- In the lower-right frame we see that QED+MSSM quantum corrections are negative at large px, and 
positive at low px, although in that case they are correcting a cross section which is already very small at 
tree-level. From the upper-right frame we learn that QED corrections are smaller than MSSM corrections, 
however they cannot be neglected at medium values for px- Leading logarithms are comparatively large, 
but in absolute terms they are smaller than in the different helicity cases. 

In the next three figures we plot differential cross sections where right handed electrons collide with left 
handed positrons. These cross sections are noticeably smaller than the case with left handed electrons. The 
reason is an accidental cancellation between Z and photon contributions in the amplitude already present at 
tree-level, as it was explained before. In Fig.QTJwe have da~^~ /dpx, which varies from 0.01 to 0.2 fb/GeV 
as pt grows from 100 GeV up to 468.8 GeV, which is the maximum value for px- QED+MSSM corrections 
vary slowly between 3%-6%, as seen in the lower-right frame. Nevertheless, from the upper-right frame we 
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dcr L ++ ' /dp T in Focus Point Scenario 
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FIG. 10: da^ + ' /dpr as a function of pt in the Focus Point scenario within mSUGRA. 
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FIG. 11: d<Tft jdpT as a function of p? in the Focus Point scenario within mSUGRA. 



see that most of these corrections come from MSSM loops and QED can be neglected. In addition, leading 
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logarithms can be seen in the left frame, and they are large and positive with the exception of very large p? 
where they can be negative. 




FIG. 12: da^ + /dp T as a function of px in the Focus Point scenario within mSUGRA. 

In Fig. [12] we have da^ + / dpr, which is a factor 1/2 smaller than the previous one. Quantum corrections 
corresponding to QED+MSSM are larger, between —3% and —8%, and interestingly enough they have 
different sign than the case with opposite chargino helicities, thus they cancel and not add to each other as 
in the case for left handed electrons. Nevertheless, these large corrections for right handed electrons act on a 
10 times smaller cross section compared with left handed electrons. If electrons are not 100% polarized (as 
expected) the pure right handed cross section will be diluted by the left handed cross section. We see from 
the upper-right frame that the QED corrections are of the same order than MSSM ones, and non-negligible. 
Leading logarithms are large, positive for small pr and negative for high px- 

Finally, in Fig. [13] we have the equal chargino helicities production cross section da^ + /dpx = 
da~^~ jdpT for the case of right handed electrons and left handed positrons. This cross section is ex- 
tremely small, never bigger than 0.02 fb/GeV. Quantum corrections are very large at small pr, but with 
cross sections so small that the chance to be observed is minimal. We include it for completness. 
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FIG. 13: da], /dpT as a function of pt in the Focus Point scenario within mSUGRA. 

VII. CONCLUSIONS 

We have calculated the complete one-loop quantum corrections to the differential cross section for 
chargino production in electron positron collisions, at 1 TeV center of mass energy, relevant for the fu- 
ture International Linear Collider. In the analysis we have included the electron and positron polarization, 
as well as the chargino helicities. 

We organize the non-QED corrections calculating the helicity amplitudes, and regularizing the ultravi- 
olet divergences with Dimensional Regularization. Our renormalization scheme includes on shell masses 
and DR running couplings. In addition, we use the dipole subtraction formalism in Bremsstrahlung con- 
tributions to cancel analytically soft and collinear divergences, with infrared divergences from virtual QED 
corrections. The well known leading logarithms remain, and we isolate them in order to better understand 
quantum corrections. 

We show numerical results for the Focus Point scenario in mSUGRA. The dipole subtraction formalism 
is designed to handle analytically the divergences, therefore the use of infrared cut-offs becomes unnec- 
essary, making the calculation numerically stable. We find important quantum corrections for differential 
cross sections with definite chargino helicities. The substantial differences in the NLO corrections depend- 
ing on the helicities of the charginos and of the initial leptons, show that the rich structure of the MSSM 
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- in particular the parity violating properties - can be investigated, firstly by polarization of the incoming 
electron/proton beam and secondly by determining (through the angular distributions of the decay prod- 



ucts) the helicities of the charginos produced j29i . 



, 40]. Since charginos are unstable, the full potential of this 
work will be fulfilled when similar radiative corrections are included for the chargino decay, and the spin 
correlations between production and decay are evaluated. Large corrections of the order of -18% and -35% 
are obtained for the largest cross sections, da^ + /dpr and da^ /dpr- They are large enough to making it 
necessary to include them in any analysis for the ILC. 
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